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Abstract 

We investigate the collective oscillations of a Bose condensate in an 
anisotropic harmonic trap. We analytically calculate the values of the trap 
deformation for which the normal modes are in resonance. In presence of 
low-order resonances, the amplitude dependence of frequencies is strongly en- 
hanced. We also study the transition to chaos of the collective modes as a 
function of the trap anisotropy. 
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The Bose-Einstein condensation (BEC) is the macroscopic occupation of the single- 
particle ground-state of a system of bosons. In the 1995 BEC has been achieved with 
dilute vapors of confined alkali-metal atoms. 1-3 The experiments with alkali- metal atoms 
generally consist of a laser cooling and confinement in an external potential (magnetic or 
magneto-optical trap) and an evaporative cooling. 1-3 Nowadays a dozen of experimental 
groups have obtained Bose condensates by using different atomic species and geometries of 
the confining trap. 

The theoretically predicted low-energy collective oscillations of the Bose condensate 4 
have been experimentally confirmed by laser imaging techniques. 5 Moreover, the dynamics 
of the collective oscillations of the condensate has been theoretically studied beyond the 
linear regime, showing strong enhancement of the amplitude dependence of frequencies in 
presence of resonances. 6 

In this paper we derive analytical formulas that give the values of trap anisotropy for 
which the normal modes of the collective oscillations are in resonance. We discuss also the 
transition to chaos of the collective modes by increasing the energy of the system, namely 
by considering oscillations of very large amplitude. Note that the chaoticity of the high- 
energy single-particle excitations of the condensate has been shown in the classical limit by 
Fliesser et al. 7 Here, we consider instead the low-energy collective excitations: they have a 
classical-like behavior but satisfy the quantum mechanical equations of a superfluid in an 
external potential. 4 

In current experiments alkali vapors are quite dilute and at zero temperature the atoms 
are practically all in the condensate. 1-3 The dynamics of the Bose condensate can be accu- 
rately described by the Gross-Pitaevskii (GP) equation 8 of mean-field approximation 



where 0(r, t) is the macroscopic wavefunction (order parameter) of the condensate, normal- 
ized to the number N of condensed atoms. The external trap is well approximated by an 
anisotropic harmonic potential 




(1) 
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U(r) = -{u lX 2 + u 2 y 2 + oj 3 z 2 ) . (2) 

and the parameter of the nonlinear term g = 4irh 2 a s /m is the scattering amplitude with a s 
the s-wave scattering length. It is useful to introduce the geometric average tUh = (ujiuj 2 uj 3 ) 1 ^ 
of the trap frequencies and the characteristic length = (h/mujh) 1 ^ 2 of the trap. In previous 
papers we have studied the properties of the Bose condensate with attractive interaction, in 
particular for 7 Li atoms. 9 ' 10 Here, instead, we consider only the case of a repulsive interaction. 

The complex wavefunction 0(r, t) of the condensate can be written in terms of a modulus 
and a phase, as follows 



0(r,t) = V / ^M)e^ t ). (3) 

The phase S fixes the velocity field v = (h/m)VS. The GP equation can hence be rewritten 
in the form of two coupled hydrodynamic equations for the density and the velocity field. 4 
If the repulsive interaction among atoms is strong enough, i.e. Na s /ah » 1, then one can 
safely neglect the kinetic pressure term in the equation for the velocity field. In this way, 
the hydrodynamic equations read 

^p + V-(vp) = 0, 

d „ f rT mv 2 \ 
m dt V+ [ U + 9P + —) =0 - ( 4 ) 

These equations are those of a superfluid in an external potential. 4 The ground-state so- 
lution (v = 0) of the hydrodynamic equations is p(r) = g^ 1 [fj, — U(r)], in the region 
where \i > U(r), and p = outside. The normalization condition on p(r) provides 
p= (hiu h /2)(l5Na s /a h ) 2 / 5 . 

An analytic class of solutions of Eq. (4) has been found found by Dalfovo et al. 6 writing 
the density in the form p(r, t) = a (t) — ai(t)x 2 — a 2 {t)y 2 — a 3 (t)z 2 in the region where p(r, t) 
is positive, and the velocity field as v(r, t) = \V\bi(t)x 2 + b2{t)y 2 + b 3 (t)z 2 }. The coefficient 
a is fixed by the normalization of the density a = (15N/87i) 2 ^(aia2a 3 Y^ 5 . By inserting 
these expressions in the hydrodynamic equations one finds 6 coupled differential equations 



for the time- dependent parameters a,i(i) and bi(t). By introducing new variables defined 
by di = muf (2(7g J 2 gig 2 g 3 )~ 1 , the hydrodynamic equations give 6j = (ji/qi and 

ft + ufq, = l - — , (5) 

Qi <?1<?2<?3 

with i — 1,2,3. It is important to observe that, using the new variables qi, the equations of 
motion do not depend on the value of the coupling constant g. In terms of the mean square 
radii of the condensate are < x\ >= (2fi/muji)q 2 and the velocities are < v 2 >= (2/t/mw,)g t 2 
(see also Ref. 6). 

The Eq. (5) are the classical equations of motion of a system with coordinates % and 
Lagrangian given by 

L = il>r 2 ?? -??)--!-. (6) 

Z i=l 9192^3 

This Lagrangian describes the collective modes of the Bose condensate for Na s /ah » 1. 
In such a case the collective dynamics of the condensate does not depend on the number of 
atoms or the scattering length. The low-energy collective excitations of the condensate are 
the small oscillations of variables q^s around the equilibrium point, i.e. the minimum of the 
effective potential X^Li qf + l/(<?i<?2<?3) ■ The normal mode frequencies Q of the condensate 
are given by the following equation 

Q 6 - 3 (ul + w\ + ufj Q 4 + 8 (u)\ uo 2 2 + u\ uol + u 2 2 ufj Q 2 - 20 u\ w% u% = . (7) 

Note that this formula has been recently obtained by using a variational approach 11 and also 
by studying hydrodynamic density fluctuations of the condensate. 12 For an axially symmetric 
trap, where lo\ = u 2 = u±, the previous equation gives 

Qi, 2 = J 2 + -A 2 ± -VlQ + 9A 4 - 16A 2 uj ± , 
y 2 2 

n 3 = V2 u ± , (8) 

where A = oo^/u^ is the anisotropy parameter of the trap. The frequencies f2 li2 are the 
m z = monopole and quadrupole modes, where m z is the third component of the angular 



momentum. Instead the frequency f2 3 is the m z = 2 quadrupole mode. 4 It is important to 
observe that the experimental results obtained with sodium vapors at MIT (A = y/8) are in 
good agreement with these theoretical values. 13 

The normal modes are functions of the trap frequencies, i.e. of the trap deformation. 
Beyond the linear regime there is a strong enhancement of nonlinear effects when frequencies 
of different modes, or of their harmonics, coincide, i.e. when the normal modes are in 
resonance. 6,14 By using a perturbative approach, Dalfovo et a/. 6 have identified 3 resonant 
values of the trap anisotropy. For such values, the numerical solution of Eq. (5) shows 
that the collective frequencies deviate from the predictions of Eq. (7) or (8) also for small 
amplitude oscillations. 

To find the BEC resonances we use a different and new strategy. In general, the resonance 
condition is given by 



riiVti + n 2 Vt 2 + n 3 fi 3 = , 



(9) 



where (ni, n 2 , n 3 ) G Z 3 \{0} is a non-zero vector of integer numbers and R = \ni\ + \n 2 \ + \n 3 \ 
is the order of the resonance. It is possible to numerically determine the values of the trap 
deformation for which the resonance condition is satisfied. Nevertheless, because in most 
experiments the confining trap has axial symmetry, 1-3 one obtains remarkable analytical 
results from Eq. (8) and (9) by imposing that one of the integer numbers rij is zero. In the 
case rt\ — 0, the resonance condition (between a m z = mode and the m z = 2 mode) gives 



A 



2(2n|-n|)n| 



(10) 



\| (5nl - 3n§)n| ' 

In the case n 2 = 0, one has the previous formula with n\ instead of n 2 . Finally, in the case 
n 3 = one finds 



5 (al + ai) _ 2 

\n 2 2 n\) 



± 



\ 



- 144 



(11) 



This formula gives the anisotropy parameter at the resonance condition between the two 



m z = modes. 
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In Table 1 are shown the values of the trap anisotropy A for low-order resonances accord- 
ing to Eq. (10) and (11). Note that the values of A for the (1,2,0) and (0, 1,2) resonances 
are that previously obtained by Dalfovo et al.. 6 As concerns resonances, the values of A 
for the experimental traps of MIT (A = 0.077) and JILA (A = y/8) are not of particular 
interest. Instead, it would be very interesting to study experimentally the system in the 
trap conditions given by Eq. (10) and (11). In fact, close to low-order resonances, even for 
relatively small amplitudes, the motion is rather complex and the dynamics can exhibit a 
chaotic-like behavior. 14 

One expects an order-chaos transition by increasing the amplitude of oscillations, thus 
by increasing the energy of the system with respect to the ground-state (minimum of the 
effective potential). 14 To check if such a scenario works in this context, we numerically study 
the chaotic dynamics of the m z = modes (preliminar results can be found in Ref. 15). We 
put qi = q 2 = q± in the BEC lagrangian. By using the adimensional time r = u±t, the BEC 
Hamiltonian of the m z = collective modes reads 

H = pi + + q 2 ± + \ql + , (12) 

where p±_ = dq±/dr and p 3 = \^ 2 dq 3 /dr are the conjugate momenta. 15 We use a symplectic 
method to numerically compute the trajectories. The conservation of energy restricts any 
trajectory of the four-dimensional phase space to a three-dimensional energy shell. At a 
particular energy, the restriction q± = 1 defines a two-dimensional surface in the phase 
space, that is called Poincare section. Each time a particular trajectory passes through the 
surface a point is plotted at the position of intersection (53,^3). We employ a first-order 
interpolation process to reduce inaccuracies due to the use of a finite step length. 

In Fig.l and Fig. 2 we plot Poincare sections of the system with A = 1.501 (corresponding 
to the (3, 5, 0) resonance) and A = y/8 (JILA experiment), respectively. In each panel there 
is a Poincare section with a fixed value of the energy of the system. At each energy value, we 
have chosen different initial conditions [q±(0),q 3 (0),pj_(0),p 3 (0)} for the dynamics. Actually 
P_l(0) has been fixed by the conservation of energy. Integration time is 400 in adimensional 
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units, that is less than 1 second (the life-time of the condensate is about 10 seconds). Chaotic 
regions on the Poincare section are characterized by a set of randomly distributed points 
and regular regions by dotted or solid curves. 

We use the adimensional parameter x, that is the relative increase of energy with respect 
to the ground-state. Fig. 1 shows that for A = 1.501 there is an order-chaos transition by 
increasing x- a chaotic sea appears at low energy (x = 2/5). On the contrary, Fig. 2 shows 
that for A = V8 at x — 3/5 the collective oscillations are inharmonic but still not chaotic. 

To estimate the chaoticity of the m z = modes, we calculate the fraction of initial 
conditions that give rise to regular trajectories. In the (X,x) plane we plot the configurations 
for which the m z = oscillations are regular. The numerical results are shown in Fig. 3. 
As previously discussed, for A = 1.501 (or A very close to such a value) the system remains 
regular only at small x; for larger x it becomes gradually chaotic. In general, the transition 
to chaos depends on x, but for A << 1 and A > 4 the dynamics appears always regular. 
Actually, for A << 1 the system is integrable because of an additional constant of motion. 16 
The MIT trap (A = 0.077) belongs to this stable region. 

We have analyzed the Poincare sections of the m z = BEC Hamiltonian for all the 
resonant values of A that are shown in Tab. 1 (left). Among the resonances, the (3,5,0) 
resonance (A = 0.888 and A = 1.501) appears the most chaotic: this is the main prediction 
of our paper. 

An important question is the following: can BEC chaotic dynamics be experimentally 
detected? The answer is positive. It is important to observe that, strictly speaking, there is 
no dynamical chaos in quantum systems with a finite number of degrees of freedom (for a 
discussion of chaos and quantum chaos in field theories see Ref. 17). Chaotic behaviour is 
possible only as a transient with lifetime tn, the so-called Heisenberg time. 18 Nevertheless, 
tjj grows exponentially with the number of degrees of freedom and consequently the transient 
chaotic dynamics of quantum states and observables can be experimentally observed in Bose 
condensates. 19 Nowadays non-destructive images of the dynamics of the condensate can be 
obtained. The radius of the condensate as a function of time can be detected and its time 
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evolution analyzed with great accuracy. Various initial conditions of the collective dynamics 
can be obtained by using laser beams or by modulating for a short period the magnetic 
fields which confine the condensate. 

In conclusion, we have studied the collective oscillations of a Bose condensate in an 
anisotropic harmonic trap. We have found the values of trap deformation for which the 
normal modes of the collective oscillations are in resonance. For such values, nonlinear 
effects, like frequency shift and irregular patterns in the time evolution of the condensate 
shape, are strongly enhanced. We have analyzed the order-chaos transition of the m z = 
collective oscillations by increasing the energy, namely the amplitude of oscillation. We have 
obtained compelling evidence that, for certain values of trap anisotropy, chaos is reached at 
very low energy. 

We hope that our calculations may stimulate new measurements on the collective oscil- 
lations of Bose condensates in anisotropic traps. 

The author thanks Prof. L. Reatto for useful suggestions. 
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Resonance Ai A 2 


Resonance A 


(1,2,0) 0.683 1.952 
(1,3,0) 0.438 3.081 
(1,4,0) 0.320 4.159 
(1,5,0) 0.255 5.226 
(2,5,0) 0.527 2.530 
(3,5,0) 0.888 1.501 


(0,1,1) 1 
(0,1,2) 1.512 
(0,1,3) 2.393 
(0,2,1) 0.454 
(0,2,3) 0.802 
(0,3,1) 0.300 



Table 1. Low-order resonances of the normal modes of the condensate and the corresponding 
values of the anisotropy parameter A of the external trap. 
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FIGURE CAPTIONS 



Figure 1. Poincare Sections of the m z = modes. From left to right: x — 1/5, 2/5 and 
3/5, respectively, x is the relative increase of energy with respect to the ground-state. Trap 
anisotropy A = 1.501. 

Figure 2. The same as in Figure 1. Trap anisotropy A = (JILA). 

Figure 3. Configurations of the (A,%) plane for which the m z = oscillations are 
regular. A is the trap anisotropy and x is the relative increase of energy with respect to the 
ground-state. Full triangles: MIT trap. Full circles: (3, 5, 0) resonance. Full squares: JILA 
trap. Open circles: all the other configurations. 
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